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o : 

_ In this work, we provide a method which aUows to compute exactly the multipoint and multi- 

■ time correlation functions of a one-dimensional stochastic model of dimer adsorption-evaporation 
, ^ ' with random (uncorrelated) initial states. In particular explicit expressions of the two-point non- 

instantaneous/instantaneous correlation functions are obtained. The long-time behavior of these 

■ expressions is discussed in details and in various physical regimes. 

\0 \ PACS number(s); 02.50.-r, 68.43.Mn, 05.50.-fq 

• One-dimensional reaction-diffusion (RD) processes of interacting particles have been extensively studied in the last 
'■^ I decade because of their relevance as examples of strongly-correlated nonequilibrium systems and their connexion with 

(]J • experimental situations [|l]-Q . 

^ I Among the RD systems, the "diffusion-limited with pair-annihilation and creation" (DPAC) model p|-p^ plays 
I . a particular role. In fact it is one of the rare nonequilibrium model for which it has been possible, in some special 
' cases, to compute some dynamical correlation functions. In addition this model carries valuable informations for 
, various experimental situations where particles diffuse and dimer can be adsorpted/evaporated [0^. Despite the 
• ' interest in the DPAC-model, not all the desirable information on the correlation functions was available so far. In 
, particular considerably less results are available for the (complete) DPAC-modcl than for the "diffusion-limited pair 
annihilation" (DPA) model (where there is no pair-creation). 

Recently, there has been a regain of interest for the study of the DPAC-model because its possible application in 
. various fields such as the experimental study of the photogrowth properties of long-lived midgap absorption band 
^1 ' in a MX chain ||^ and in interdisciplinary studies In particular, it has been shown that the autocorrelation 
functions of the DPAC-model provide valuable informations on the relaxation of biological dimer adsorption |^ . In 
this work we consider the (free-fermion) DPAC and DPA-model and obtain results which remained inaccessible so far. 
In particular, we explicitly compute the exact and complete expression of the noninstantaneous two-point correlation 
functions for random initial conditions and then analyse the long-time behavior of the latter. 

We consider a periodic lattice of L sites (without restriction, L is assumed to be even) on which an even number 
of (classical) particles interact. Each site is either empty or occupied by a particle at most (because of the hard-core 
interaction). When a particle and a vacancy are adjacent to each other, the particle can jump to the right with 
, a rate h' or to the left with a rate h. When two particles are adjacent, they can annihilate in pairs with a rate 
' e. In addition, when two vacancy are adjacent, a pair of particles can be created with rate e'. We now adopt the 
' so-called stochastic Hamiltonian formalism (see e.g. [Q and references therein). To do this, at each of the lattice 

• sites, we associate to a particle (vacancy) a spin-i down (up) . In so doing the master equation governing the 
dynamics of the model can formally be rewritten as an imaginary-time Schrodinger equation for a quantum spin- 
chain: [d/dt)\P{t)) = —H\P{t)), where \P{t)) = X]{n} ^({"li OK"-}) describes the state of the system at time t 
(the sum runs over all the 2^ configurations {n}). Performing a standard Jordan- Wigner transformation and then 

'"^ a Fourier transformation, the stochastic Hamiltonian can be recasted in a fermionic representation p|-p^]. We also 
', define the "left vacuum" {x\ = J2{n}ii^}\- '^^'^ probability conservation yields (xl-^ = 0. 
^ ■ Exact solution of the DPAC-model is only possible in the free-fermion case ||^-[l^ , and therefore, with 7 = 
e-\- e' — {h-\- h'), one has to impose: 7 = 0. 

Hereafter, we always consider that the constraint 7 = is fulfilled. In this case, the stochastic Hamiltonian reads: 



■ H = '^\uj{q)a\aq+ u)*{q)a^_^a-q + 2sixiq(^eaqa-q + e'a^_qa^^ + e' L (1) 

where and Oq are usual fermion operators. In addition, uj{q) = c — bcosq -\- ivsinq, with 6 = e + e', c = e — e', 
v = h'-handq = ±7r(2Z - 1)/L, Z = 1, 2, . . . , L/2. 

We consider translationally invariant and uncorrelated random initial states |po) with an even number N of particles, 

of density po — N/L, i.e. \pq) — i p^^ ) ' fermion reformulation, an important property of the left 

vacuum follows (in the subspace with an even number of particles that we consider here): (x|aj = cot(|)(^|a_q. 
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The free-fermion character of the stochastic Hamiltonian allows the computation of the following zero-time 
correlators ||9|, with n = po/{l — po): 



^2 cot {q/2) 
1 + /i2 cot2 (g/2) ' 

(a^jOg^ . • •ag2„-ia<Z2„>(0) = ^ 5(7r)(a,^(^,a5^(,))(0) . . . (a?,(2„-i)a9,(2„))(0), (2) 



where the sum is over all the permutations tt of the indices {qi, q2, ■ ■ ■ , q2n}, with the constraints vr(l) < 
7r(2); . . . ; 7r(2n — 1) < 7r(2n). Each permutation tt has a signature 5(7r) . 

It is advantageous at this point to introduce the "pseudo-fermion" operators |^ = cos 9qaq + a sin Oqa}_^ and 
^+ acos^qfll^ + sm6qa-q , with = —9q. Although they are not adjoint each other Q, these operators fulfill 
the canonical anticommutation relation {S,q,^q,} — Sg g/. The probability conservation yields (xlCq'^ = 

It has been shown lloll that if one chooses the following parameters tan(26'„) = "^^^^^tS ■ o? — J t' It. the stochastic 

free-fermion Hamiltonian (|^) is diagonal in the pseudo-fermion representation and reads: H = Aq^^^^, with Ag = 
b — ccosq + TO sing. Because of this diagonal representation, the pseudo-fermion operators evolve according to 
Ut) = e-^'*Cg(0) and = e^«*^+(0). 

Using the fact that, the pseudo-fermion operators are linear combination of fermion operators Oq and aj, and 
using the expression (Q) as well as the property of the left vacuum, it follows that: 

i^q^q' > (0) = ^"^^gcos^?' ^^^^^^ ^ ^ ^2 gjj^ (a%a^_g, ) (0) + sin 9q, cos 0, (agal^, ) (0) + sin 9q cos {al^Oq, ) (0) (3) 



/(g) ^ (^,e-,)(0) = - [1 + 2^.,^ + 2/^.2] , where i^q ^ ^ecot (g/2), (4) 



For the sequel it is useful to compute {^q£,q'){0) — {x\^q^-q\po)Sq' ,-q- Therefore we introduce the following function : 

fi'^evq cos^{q/2) 

One can also check that 

where we adopted the same notations as in (|^). 

Let us now sketch a four-step procedure which allows to compute explicitly the multipoint and multi-time correlation 
functions {n.j^{ti) . . .nj^_,{t^-i)nj^{t^)) = {x\nj,e-"''*^-*^^ . . . nj^_,e-"'-*'"-^-*'^'^nj^e-"^'"\po): 

(i) One first has to write the expression of the correlation functions in the Fourier space: 

(ii) One has then to rewrite the expression in the pseudo-fermion language. This achieved with help of 

ajflg/ = cos cos 6*5/^+^5/ - cos6'gSin6lg/^+^+^, - sin 61, cos 6'q'^_g^_g' + sin6'gSin6'g/^_g^+^,. (7) 

(iii) Using the fact that H is diagonal in the (pseudo-fermion) representation [with (|^)], one extracts the time- 
dependence of terms appearing in the pseudo-fermion rewritten expression (^). As an example, we have: 

(iv) Finally, the zero-time correlation functions of pseudo-fermion operators appearing in the r.h.s. of (|^), after the 
steps (i)-(iii), are computed with help of the Wick factorization (||) and using (x|C^ — 0. 

This general four-step procedure provides a systematic method to obtain explicitly, starting from homogeneous 
random initial conditions, the multipoint and multitime correlation functions of the free-fermion model under con- 
sideration here. It is however important to notice that the computation of each of these quantities leads to rather 
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complicated technical difficulties. Let us mention that using the domain-wall duality and with help of the generating 
function studied in , we can compute the stationary multipoint correlation functions of the DPAC-model from 
the spin-spin correlation functions of the one-dimensional Ising model with a generalized (biased) Glauber's dynamics 
P, p^p^ . In fact we can show > . . . > ji) that {rij^ . . . nj^J(oo) = [p(oo)]™ where p(oo) = /(-/e -I- is 
the stationary density of particles. In addition, for an homogeneous system with initial density po = 1/2 of particles, 
because of the quadratic form of the generating function, the expressions of spin-spin correlation functions of the dual 
of the DPAC-model are Pfaffians In this case it is therefore possible, via the domain-wall duality to 

sort out the technical complications and explicitly compute the instantaneous multipoint correlation functions. As 
an example for the DPA- model [e ~ h + h' , e' ~ 0)^ the long-time behavior {bt ^ 1) of the three-point correlation 

functions read 8(n,(t)n,+,, (t)) « ~;3"'^'^"'^^ ^ith > n |§. 

Because of the general technical problems inherent to the computation of the correlation functions, the above- 
mentioned systematic four-step procedure is in particular useful to take into account random initial conditions, which 
affect the long-time dynamics of the non-universal relaxation (when all the rates h, h' , e and e' are > 0, see below). 

To illustrate the difficulties which appear in computing the multipoint and multitime correlation functions (from 
uncorrelated but random initial states), as well as their importance, one can point out the work of Derrida and Zeitak 
[ prt , where these authors obtained the universal distribution of domain sizes of one-dimensional Potts model with zero 
temperature Glauber dynamics. This was achieved using the properties of coalescing random walkers to compute the 
probability of having the same value at time t for N spins located at N distinct and ordered sites, which is related 
to the distribution of domain sizes. The authors also studied the domain-walls dynamics and thus considered the 
following (free-fermion) RD model ^ 0A, A A A and A A — ^ 0, with reaction-rates 1, (q — 2)/(q — 1) and 
l/(q— 1), respectively. For this RD model, with random (but uncorrelated) initial conditions, the authors of 
computed the density of particles and the instantaneous two-point correlation functions. It has to be noticed that 
there exists a similarity transformation (see e.g. fe| and references therein) which maps the DPA-model studied here 
(with e' ~ 0) onto the RD model considered in [p7|. 

Another important problem where the (stationary) multipoint correlation functions of a free-fermion model play 
a relevant role is the computation, for the q-state Potts model in the zero-temperature Glauber dynamics, of the 
exact persistence exponent which gives the fraction of spins pL{<i) which have never flipped [T^Jl9|] . To compute 
Pl(q.)i the authors mapped the problem onto an exactly solvable (free-fermion) RD model: A0 «-> 0A (with reaction- 
rate 1) and AA — > A (with reaction-rate 2). In addition it is assumed that a "source" ensures that the origin of 
the (periodic) lattice is always occupied. Starting from an uncorrelated but random initial state (with the initial 
site always occupied) denoted |P"(0)) , the problem of finding PL(q) reduces to the computation of the following 
multipoint correlator p8| : pL(q) — hmt^oo(0|(l + ul) ... (1 + ai)e^^ *|P"(0)), where the a^-'s are fermion operators, 
|0) is the vacuum (aj|0) = 0) and H' is the stochastic Hamiltonian associated to RD model considered in ||l^ , pj|| . 

To be specific we now focus on the computation of the connected non-instantaneous two-point correlation functions 
of the DPAC-model for random initial conditions |po)- 

Following the above-mentioned four-step procedure [(i)-(iv)], adopting the notation (j)q = Rc(Ag) — b ~ ccosq, we 
obtain, in the thermodynamic limit [L — )■ oo): 



Cr{tM) = {nj+r{t + to)nj{to)) - p{t + to)p{to) = ee ( f — sin[9r - wising] ^e""^'* 

. > f r dq , . ,sin^(q/2) ^ t\ f T dq , . ^cos^(q/2) 
-h4ee' / — cos [gr - wising] llL^e""^'* / — cos [gr - wising] yiLJ.Q''l',t 

+ 2e'V7^( r ^sm[qr-vtsmqf^e~^A ( f ^ sin[gr - wi sin g] ^^^^M^/(g)e-^«(*+2*'')^ 

V^O 7^ <Pq J \Jo I" 4>q J 

r ^cos[qr^vtsmq]^^^-^^e-^A ( T ^ cos[gr - wi sin g] ^/(g)e-'^'(*+2*°)') 

\Jo 4>q J \Jo 0, y 

-2eV7^( r ^sm[qr-vt sin qf-^e-'^ A ( T ^ sin[gr - w<sing]^^^^l^/(g)e-^'(*+2to) 

\Jo (Pq J \Jo 



q 

;„2 



4«'f r^?«sin[gr-wising]^^^i^/(g)e-^«(*+2*'')^ ( C ^ sm[qr ~ vt sin q f" '"'"' fiq)e 

\Jq (Pq J \Jo ^ 'Pq 

ee'f r^^cos[gr-wising]^/(g)e-^'(*+2*°)) (9) 



09 



3 



In this expression p{t) ~ ^/^^^[/^ ~ lo ^^^f^l)^ designates the (translationally-invariant) density of 
particles at time t. The latter has been previously studied (for the initial states |po)) in 0]- 

The equation is the main result of this work and provides the complete expression of the noninstantaneous 
two-point correlation functions of the (free-fermion version of the) DPAC-model. From the latter, it is clear that one 
can also obtain the instantaneous two-point connected correlation functions Crit = 0,to > 0)i and one can check , as 
expected from the general properties of the DPAC-model that the latter do not depend on the bias v. In addition, it 
is clear that the expression (|^) also includes the instantaneous and noninstantaneous connected two-point correlation 
functions of the DPA-model, where e' = and 6 = c > 0. It is worthwhile to notice that the noninstantaneous 
correlation functions (^) [with t > and v ^ {)] depend on the sign of r : Cr{t,tQ;v) = C-r{t,to;~v). Conversely, 
the instantaneous correlation functions [with t = and to > in (H) ] do not depend on the sign of r. Let us also 
stress the fact that setting in (^) /i = oo (i.e. po = 1), fj, = (i.e. po — 0) or p = 1 (i.e. po = 1/2), we recover results 
obtained in |-|l||. 

To proceed with a long-time study of the expression (^) , we have carried out a systematic asymptotic expansion 
of the integrals appearing in (^, in which the small q regime of integration dominates. Hereafter, we analyse two 
different regimes and distinguish the case with pair creation (i.e. with e' > and c 7^ 0) from the case with only 
(asymmetric) diffusion and pair-annihilation (i.e. with e' = and e > 0). The cases where c — (i.e. e = e', /i = /i'; 
and e = e', h + h' = 2e) have already been studied in [^. 

l.i) We first consider the case where e' > and e > in the regime where et, eto ^ 1 and e't, e'to ^ 1 [with c > 0]. 
In this situation the main contribution to the noninstantaneous correlation function arises from the second term of 
the r.h.s. of (|). We thus obtain {v ^0,c> 0): 

C(t,to)^ u=ir-vt)yct (10) 

It is clear from ( |l0| ) that in this regime the late behavior of the noninstantaneous correlation functions Cr{t,to) only 
depends on the time t (and not on to). We notice the non-trivial effect of the bias v ^ through the parameter u. In the 

e exp[— 4e'i] 



absence of the bias and for r < 00 (i.e. for the autocorrelation functions) we obtain: Cr{t,tQ;v — 0) ^ ~ ^ 

In this regime we now focus of the long-time behavior of the instantaneous correlation functions Cr{t) (obtained 
setting t = in (|^) and relabelling as the variable t) . The cases po = I and pQ = having been studied previously 
[Cr{t; Po = 0, 1) cx e"^*^ *t~'^ , 1^' = 3/2 for po = 1 and i^' = 1/2 for po — here we focus on the case of random 

initial states, i.e. with < p. < 00 [and po 7^ pioo) = \fe! /(Ve' + v^), which would correspond to the trivial case where 
f(q) = 0]. The main contribution to Cr{t) comes from the third and fourth term of the r.h.s. of (^. Introducing the 
parameters: 




(11) 



Tre ^^'*J^{p,T,e,e') 



where the rather complicated expression of the amplitude J^{p, r, e, e') reads: 

:F{h, r, e, e') = {2p\ee' f' - 2e''V77{pr)' + ee' [^e + p^[2r^^e- SV^]) } 



{6/i'(6e')'/' - /i'(2r2 + l)e' W + ee' (z^^e + [(2^2 + l)^^e- 9V^] ) } (13) 



12ee'^p 

It is remarkable that conversely to the instantaneous correlation functions ([T2[), which amplitude (^|) depends 
on the initial state through the parameter < /i < cxd, the long-time behavior of the noninstantaneous correlation 
functions ( p^ do not depend on po- This is due to the fact that the second term of (||) does not depend on f{q). 

l.ii) Another interesting asymptotic regime to investigate is the one first studied by Torney and McConnell [psj, 
where one considers initially very diluted systems, i.e. po ~ p ^^ but keeps the products ep^t, epQto, e'pgt, e'p^o 
fixed and finite, with et, eto 3> 1 and e't, e'to 3> 1 [and c > 0]. 
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In this regime, the noninstantaneous two-point correlation functions Cr(t, tp) are stih dominated by the second term 
of (^) and thus the asymptotic {v ^ 0) decay of Cr{t,to) is still given by 

The situation is however different for the instantaneous correlation functions [because the third and fourth term of 
depend on /((?)]. With help of (|Tl|), we obtain (r < oo): 



Crit) 



p3e-4^'*e(po,r,e,e') 



e4po^tErfc(2poVa) 



(14) 



_2(p27ra)i/2 

where Erfc(z) denotes the usual complementary error function. The amplitude G{pq, r, e, e') has the following form: 



g{po,r, e, e') = - Cq) [pq ' - + e/e'] - Sqv^ r/pl - 10r{2r^ + 1) - ^ 



(15) 



We now pass to the case of the DPA-model, where e' = and b = c = e > 0. Again, we distinguish two different 
regimes. 

2.i) We first consider the regime where bt ^ 1 and btp 3> 1. In this situation, the main contribution to the long-time 
dynamics arises from the fourth and the last terms of (|) and we obtain : 



Cr{t, to;e' = 0) 



exp 



(r~vty to+t 
bt 2to+t 



2Trby/t{t + 2to) 



t 



t + 2to 



2to+t 



(16) 



It has to be noticed that, according to (|16[), when r > and v > (or, r < and v < 0), Cr{t, to) has a local maximum 



too 



(a "peak") at time tp = r/v. When bt '3> bto '3> 1, we recover the result ||7|,|ll|: Cr{t, to; e' = 0) S^m^" 

In this regime, the main contribution to the instantaneous correlation function arises from the last term of (|^), and 
we obtain the following result: 



Crit;e' = 0) w -l/Airbt, 



(17) 



where the minus sign manifests the fact that the long-time dynamics of the DPA-model is dominated by anticorrelation, 
due to the pair-annihilation of the particles. The results (|ig ) and ( p^ and the fact that the latter do not depend on 
po confirm, for random initial case, the universal character of the DPA-model in this regime. 

2ii) We now consider the low-density regime of the DPA-model, where po ~ M ^ 1 and et, eto ^ 1) with ep^t and 
epgto finite. Also in this regime the main contribution to Cr{t, tQ;e' = 0) arises from the fourth and the last terms of 



and one has the long-time behavior (w = and r < oo) 



Crit,to;e' = 0) « poe2''o''(*+2*«)Erfc {2poVKto + t/2) 



1 



^^g2p„6(t+2to)E,fc i^2po^/b{to + t/2) 



y/2nbt 

For the instantaneous correlation functions, we obtain the following long-time behavior: 

2 



Cr{t;e' = Q) 



'Po 



Erfc(2poV6t) 



(18) 



(19) 



where the anticorrelated character of the DPA-model clearly appears. 

Despite the fact that the parameter c = e — e' can take negative values, so far we have always considered the case 
where c > 0. With help of the similarity transformation B — Ilj^i '''J j where crj is the usual Pauli's matrix acting on 
the site j, we show that the case where c < is directly related to the one where c > 0. In fact, according to B, the 
(free-fermion) DPAC stochastic Hamiltonian H = H{b,c,v) is mapped onto BH{b,c,v)B^^ = H{b,—c,—v) and the 



initial state \po) is mapped onto B\po) = |1 ^ Po)- Therefore, we have: Cr{t,to), 



Cr{t, to) 



h.c.v\ 1 — Po ' 



In summary, in this work we sketch a four-step procedure which allows the explicit computation of the multipoint and 
multitime correlation functions of the free-fermion DPAC-model starting from random (uncorrelated) initial states. 
We then specifically compute the noninstantaneous/instantaneous two-point correlation functions in the presence as 
well as in the absence of the pair-creation term. When all the reaction-rates are positive, the dynamics turns out to 
be non-universal and the long-time relaxation is exponential (with a subdominant a power-law factor): the amplitude 
of the instantaneous two-point correlation functions depends on the initial density po and is explictly determined. 
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In the absence of the pair creation, i.e. when e' = and h + h' = e > 0, the dynamics turns out to be universal 
(in the regime where po is finite and bt ^ 1, bto 3> 1) and there is a power-law relaxation. The effect of the bias 
V = h' — h 0, only appears in the noninstantaneous correlation functions and can be absorbed (for e' > 0) in a 
Galilean transformation, as noticed in |7||Tll] in considering the DPA- model (in these previous works, e' = 0). 

To close this work, it is natural to wonder what is the effect on the dynamics of the restriction 7 = 0. In fact, it is 
by now well established on the basis of numerous consistent numerical results [^||^ Jll| , [l^ , and from comparison with 
experiments ||^,|^ , that the results obtained for the free-fermion version of the DPAC-model give a qualitative picture 
which is still valid when 7 7^ 0. One can therefore expect that the results obtained in this work could have a general 
validity and, in particular, a direct relevance for recent interdisciplinary studies Q|. 

We are grateful to P.-A. Bares, M. Michalakis and L. Zuppiroli for discussions. The support of the Swiss National 
Foundation is acknowledged. 
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